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Abstract 

In many models of the Calvin cycle of photosynthesis it is observed 
that there are solutions where concentrations of key substances belonging 
to the cycle tend to zero at late times, a phenomenon known as overload 
breakdown. In this paper we prove theorems about the existence and 
non-existence of solutions of this type and obtain information on which 
concentrations tend to zero when overload breakdown occurs. As a start¬ 
ing point we take a model of Pettersson and Ryde-Pettersson which seems 
to be prone to overload breakdown and a modification of it due to Pool- 
man which was intended to avoid this effect. 


1 Introduction 

Photosynthesis is one of the most important processes in biology and a variety of 
mathematical models have been set up in order to describe it. These have mainly 
been concerned with the part of photosynthesis known as the dark reactions, 
also known as carbon fixation or the Calvin cycle. In the simplest pictures 
of this process (see e.g. DP) the substances included in the description are 
the five carbohydrate phosphates RuBP, PGA, DPGA, GAP and Ru5P. In 
the properties of some models of this type were considered, together with a 
model which in addition includes the concentration of ATP as a variable. In 
[T^ dynamical properties of solutions of these models and related ones were 
studied. It was found that for many of the models there are solutions where the 
concentrations become unboundedly large at late times although they remain 
finite on all finite time intervals (runaway solutions). The one exception is 
the model including ATP, for which it was shown that all solutions are globally 
bounded. At the same time it was shown for all the models considered that there 
are large classes of initial data for which the corresponding solutions are such 
that the concentrations of all carbohydrate phosphates tend to zero as time tends 
to infinity. In other words all these concentrations become arbitrarily small at 
late times. In [T^] no biological interpretation was offered for this behaviour. 
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Both these phenomena, where the concentrations become arbitrarily large or 
small, might be taken as indications that these models are inappropriate. It 
suggests that it would be worthwhile to examine alternative models to see if 
they exhibit similar behaviour. 

In this paper we study some models of the Calvin cycle which incorporate 
more aspects of the biology. In the simple models the five substances included 
are linked by reactions forming a cycle. In the models considered in what follows 
this simple circular topology of the reaction network is replaced by a more com¬ 
plicated branched one. The starting point is a model introduced by Pettersson 
and Ryde-Pettersson in [5]. We call this the Pettersson model. The unknowns 
are concentrations of substances in the chloroplast, where the Calvin cycle takes 
place. Most of the reactions in the model convert some of these substances into 
others. It also includes some transport processes, where substances are exported 
from the chloroplast to the surrounding cytosol and a process in which they are 
stored as starch in the chloroplast. It was found in that for certain parameter 
values, corresponding to high concentrations of inorganic phosphate in the cy¬ 
tosol, solutions of this model can undergo a process called overload breakdown 
where the export of sugars cannot be maintained. Since this behaviour may 
be biologically unrealistic Poolman m introduced a new model, which we call 
the Poolman model, in an attempt to avoid solutions of this type. Overload 
breakdown corresponds to a situation where more sugar phosphates are being 
exported from the system than can be produced and to combat this Poolman 
introduced an extra reaction describing starch degradation, i.e. the production 
of sugar phosphates from starch. This is the only difference between the re¬ 
action networks used in the Pettersson and Poolman models. There is also a 
difference in the kinetics which describe the mechanisms of certain reactions. 
Yet another choice of kinetics, which has the advantage of simplicity, is mass 
action kinetics. The models obtained from the two networks by applying mass 
action kinetics will be called the Pettersson-MA and Poolman-MA models. For 
all these models it is easy to see that all solutions remain bounded, due to the 
conservation of the total amount of phosphate. They have no runaway solu¬ 
tions. In what follows we investigate under what circumstances these models 
with mass action kinetics admit solutions where some concentrations tend to 
zero at late times and to what extent these conclusions can be transferred to 
the original Pettersson and Poolman models. 

The paper proceeds as follows. In Sect. 2 the various models are introduced 
and the relations between them are described. The question, which combina¬ 
tions of the variables in the system might tend to zero at late times is examined 
in Sect. 3. In other words, necessary conditions for a point on the boundary 
of the state space to be an w-limit point of a positive solution are obtained. 
It is shown in Sect. 4 that in the case of the Pettersson-MA model there are 
large classes of solutions for which many concentrations converge to zero as 
t ^ oo and it is investigated to what extent this is prevented by moving to 
the Poolman-MA model. The set of substances whose concentrations may tend 
to zero gives a new picture of the details of overload breakdown. The main re¬ 
sults on the Pettersson-MA and Poolman-MA models are contained in Theorem 
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1. Sect. 5 is concerned with generalizing some of these results to the original 
Poolman model. Conclusions and an outlook are given in the last section. 

This paper is based in part on the master’s thesis of the first author [S]. 


2 Basic definitions 

This section gives the definitions of four models of the Calvin cycle mentioned 
in the introduction, the Pettersson-MA, Poolman-MA, Poolman and Pettersson 
models. The underlying reaction network is given in Figure 1. 



Figure 1: diagram of the network 


In the case of the Pettersson model k ^2 = 0 and one of the reactions can be 
dropped. In the case of the Poolman model /C 32 7 ^ 0. The system of ordinary 
differential equations obtained by applying mass action kinetics in all reactions 
will now be presented. To our knowledge this model has not previously ap- 
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peared explicitly in the literature although it does occur implicitly in [1]. In 
that reference the author investigates the applicability of some theorems of 
chemical reaction network theory to this system and since the hypotheses of 
these theorems include mass action kinetics the only logical interpretation is 
that he applied these techniques to the mass action system arising from the 
given network. The unknowns in the system are the concentrations of RuBP, 
PGA, DPGA, ATP, GAP, P;, DHAP, FBP, F6P, E4P, X5P, SBP, S7P, R5P, 
Ru5P, G6P and GIP. The equations are 


dXRuBP 

dt 

dxpGA 
dt 

— k^XpGAXATP — k 28 XpGA, 
dxDPGA 


— —kiXpuBP + k22Xpu5PXATP^ 

= 2kiXpuBP + k2XDPGA{cA — XATp) 


dt 


= -k2XDPGA{cA — XATp) + k^XpGAXATP 


+ k4XGApXPi — k^XppGA, 
dXATP 


( 1 ) 


( 2 ) 


(3) 


dt 


= k2XDPGA{cA — XATp) “ k^XpGAXATP - k22XRu5PXATP 


— k27XGlpXATP + k3iXp^{cA — XaTp) 
dxGAP 


(4) 


dt 


= -kiXGApXPi + k^XDPGA + k^XpHAP — k^XGAP + k^XpBP 


— kgXGAPXDHAP + knXpAPXXbP “ ki2XFepXGAP + kiQXx5PXR5P 

— knXSTPXGAP — k29XGAP, (5) 

dxp. 


dt 


= —kiXGApXPi + kgXpPGA + kiQXpBP + kigXSBP 


+ 2k27XGlpXATP + k28XpGA + kggXGAP + ksgXDHAP 
-k3iXp,{cA - XATp) - k32Xp,, 
dxoHAP 


dt 


= —kQXDHAP + k7XGAP + k8XFBP — kgXGAPXpHAP 


+ ki3XSBP — kiAXDHAPXEAP “ k3oXDHAP 
dxFBP 


dt 

dxFGP 

dt 


= -ksXFBP + kgXGAPXpHAP - kigXFBP, 
= kigXFBP + kiiXEAPXXbP — kigXFGPXGAP 


+ k23XG6P — kgAXFGP, 
dxEAP 


dt 


= —kiiXEAPXXbP + kigXFQPXGAP 


+ ki3XSBP — kiAXDHAPXEAP, 

dxx5P 


dt 


= —kiiXEAPXXbP + ki2XF6PXGAP — kigXX5PXR3P 


+ ki7XS7PXGAP + kggXRubP — kgiXxbPi 
- = —ki3XSBP + kiAXDHAPXEAP — ki^XSBP, 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


dt 
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dxs7P 

dt 

dxRzp 

dt 

dxpubP 

dt 


= ki5XSBP + kiGXX5PXR5P — knXSTPXQAP, 

= —kiGXX5PXR5P + ki7XS7PXGAP + kisXRubP — kigXRGP, 
= —kisXRuGP + kigXRGP — kggXRubP 


+ k2lXXbP — k22XRu5PXATP 

dxGSP 


dt 

dxGlP 

dt 


— —k2zXG6P + k2AXFGP + kg^XGlP — kgeXGGP, 

= —k25XGlP + k2eXG6P — k27XGlPXATP + k^gXp^ 


(13) 

(14) 

(15) 

(16) 

(17) 


where xx denotes the concentration of any substance X and the ki are constants, 
the reaction constants. These are all positive with the possible exception of k 32 - 
When it is zero the equations define the Pettersson-MA model and when it is 
positive they define the Poolman-MA model, ca = xadp + xatp is the total 
concentration of adenosine phosphates. It is constant in time and this fact has 
been used to eliminate xadp from the evolution equations. The concentrations 
of substances in the cytosol in the diagram of the network, which carry the 
subscript ‘ext’, are not included in the model. Their concentrations have been 
assumed fixed so that they are not dynamical variables. The same is true for 
C02, NADP, NADPH and starch. The solutions of biological interest are those 
which are positive (i.e. all concentrations are positive). Solutions which are non¬ 
negative but not positive are of interest as limits of the biologically applicable 
ones. Solutions which start positive remain positive and those which start non¬ 
negative remain non-negative. This is a consequence of Lemma 1 in the next 
section. Let S denote the positive orthant in the space of concentrations and S 
its closure. Then S and S are invariant under the evolution. The total quantity 
of phosphate in the system is 


Cp = 2xruBP + XPGA + 2.XDPGA + XgAP + Xp^ + XpHAP 

+2xfbP + XF6P + Xe4P + XX5P + 2.XSBP + XS7P 

+XRGP + XRu5P + XG6P + XGIP + SxATP + 2.XADP- (18) 

It is conserved for biological reasons and of course it follows directly from the 
evolution equations that the time derivative of this quantity is zero. In par¬ 
ticular it is a bounded function of time. Since all substances occurring in the 
system contain phosphate it follows that all concentrations are bounded and 
that solutions of the system exist globally in time. 

In the original Pettersson model a distinction is made between fast reversible 
and slow irreversible reactions when specifying the kinetics. The irreversible 
reactions are RuBP -A PGA, FBP -A F6P, SBP -A S7P, Ru5P -A RuBP, 
ADP + P^^ ATP, GIP + ATP -a ADP + P„ PGA -a Pi, GAP -a P^ and 
DPGA -A Pi. In the original Poolman model these reactions are given the same 
kinetics as in the Pettersson model. These kinetics are chosen on the basis of 
experimental data. In the Poolman model the remaining reactions are given 
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mass action kinetics. The evolution equations can be expressed in terms of 
reaction rates Vi without the kinetics being fixed. The result is 


dxpuBP 

dt 

(19) 

dxpoA „ 

= 2vi V 2 VPGA, 
dt 

(20) 

dxDPGA 
dt = 

(21) 

dXATP 

=Vie V2 Vl3 Vst, 

dt 

(22) 

dxGAP 

=V3 V4 V7 Vio VQAP, 

dt 

(23) 

dxpHAP 

= V 4 Vs VDHAP , 

dt 

(24) 

dxpBP 

dt 

(25) 

dxpQp 

= ve vr vi4. 

(26) 

dxpAP 

dt 

(27) 

dxxbP 

=V7+Vio V12, 

dt 

(28) 

dxsBP 

dt 

(29) 

dxs7P 

dt 

(30) 

dxpzp 

dt 

(31) 

dxpu5P 

=Vii+Vi2 Vi3, 

dt 

(32) 

dxG6P 

dt 

(33) 

dxGlP 

= Vi5 Vst + ni7 

dt 

(34) 


—^ = -yg + i;g + r;g + VPCA + VGAP + VDHAP + 2nst “ ^16 “ ^'l7■ (35) 
at 


These equations are identical to those in except for the fact that a rate vn has 
been added to accommodate the degradation of starch in the Poolman model 
and that the equation for the concentration of inorganic phosphate has been 
included explicitly. In this last equation was omitted since it can be computed 
from the other concentrations using the conservation law for the total amount of 
phosphate. In [S] the slow reactions correspond to the Vi with i = 1,6, 9,13,16 
and the rates of the export reactions upcA, 'Cgap and udhap- Poolman also 
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takes the additional reaction he introduces to be a slow reaction, with rate m 


vir 


_ Virxpj _ 

XPi + + K~^Ixgip) 


(36) 


The expressions for the Vi in the Pettersson-MA and Poolman-MA models can 
be read off by comparing these equations with the evolution equations for those 
models given earlier. This gives the expressions for the Vi in the Poohnan model 
in the case of the fast reactions. The rates of the slow reactions in the Poolman 
model apart from vn are taken from those in the Pettersson model and are as 
follows 


Here 


ViXpuBP 

Vi = -, 

xruBP + Ml 

(37) 

VqXfBP 

"^6 — 

xfbp + Kmei^ + + K~q\xp.) ’ 

(38) 

VexsBP 

^9 - /-I T^—1 \ ’ 

XSBP + Km9{^ + Xpj 

(39) 

Vi3XRu5PXATP 
''''' “ Mi3 

(40) 

_ yieXADpXPi 

{xadp + Armi6i)(a^Pi + KmiQ2) ’ 

(41) 

VexXpGA 

= NKpc^ ■ 

(42) 

VexXGAP 
^gap - . 

(43) 

VexXoHAP 
- NKdhap ’ 

(44) 

VstXGlpXATP 

Vaf = - . 

{XGIP + Kmstl)Mst 

(45) 


- ^iilb^ADp)], 

Kmst2^Pi 


\ -“-ilJ. -“-ilC 

= [xru5P + 

X [a;ATp(l + K~^\^xadp) + ^1^132(1 + K~^\^x 

Mst = (1 + K^J^Xadp){xaTP + Armst2(l + Ktt) 

'■{KastlXpGA + Kast2XFeP + KastsX FBp)~^), 

" — 1 4- n 4- -1 \ f 4- 4_ ^GAP XdhAP ^ 

1 + Kdhap) ■ 


X 


N 


(48) 

(49) 


It been pointed out in [2] that the expression for Vst in [9] is incorrect and 
in the expression for M^t given above we have used the replacement proposed 
in P]. In the original Pettersson model the fast reactions are assumed to be 
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at equilibrium which leads to a system of differential-algebraic equations. The 
latter model will not be treated further in the present paper. 

In the Calvin cycle most reactions conserve the number of carbon atoms. 
There are, however, some inflow and outflow reactions for which the number 
of carbon atoms which are in the chloroplast and not stored as starch is not 
conserved. Hence the total number of carbon atoms in the substances included 
in the model is not conserved. It is nevertheless useful to consider the following 
small modification of the total number of carbon atoms. Define a quantity Li 


by 


5Li = 5xruBP + -j^XpOA + 3 xdPGA + ^XqAP + ^XdHAP + GxpBP 
+ 6xFeP + 4:XeAP + 5xx5P + 7xsbP + 7xs7P + 5xr5P + 5xru5P 
-\-QxgqP + &XG1P- 


(50) 


This is inspired by a Lyapunov function constructed by trial and error in [12] 
which is related in a similar way to the total number of carbon atoms. Its time 
derivative is given by 



(51) 


Gk27XATPXGlP — ?>k29XGAP “ ‘^k^^XpHAP + k 32 Xp. 


for the Pettersson-MA and Poolman-MA models. 

In the Pettersson-MA model if it is assumed that the parameters are such 
that ksCA < 5^28 then Li is a Lyapunov function and it is strictly decreasing 
for positive solutions. It follows that in this case all w-limit points of a positive 
solution are on the boundary of S and that, in particular, there are no posi¬ 
tive stationary solutions. This parameter restriction is the direct analogue of 
one found to play a role in the dynamics of the simple model including ATP 
considered in Section 6 of m- 

3 Potential w-limit points 

In this section information will be obtained on the location of w-limit points 
of positive solutions of the Pettersson-MA and Poolman-MA systems. Many of 
the arguments use the following simple lemma. 

Lemma 1 Consider an ordinary differential equation of the form u{t) = —a{t)u{t)+ 
b{t) where a and b are non-negative continuous functions and a solution u{t) 
which satisfies u{to) > 0 for some to- Then if b{ti) > 0 for some ti > to it 
follows that u{ti) > 0. 

Proof Suppose first that u(to) > 0. Then u remains positive for t slightly larger 
than to and as long as it does so the equation can be rewritten as ^ (log u) = 

—a + Thus (log It) > —a and logu cannot tend to —oo in finite time. 

Hence u cannot tend to zero in finite time and u remains strictly positive. It 
then follows using the continuous dependence of solutions on initial data that 
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if u starts non-negative it stays non-negative. Now suppose there were a time 
ti > to with u{ti) = 0 and b{ti) > 0. Then u{ti) > 0. This implies that u{t) 
is negative for t slightly less than ti, in contradiction to what has already been 
proved. This completes the proof of the lemma. 

It is also often useful to apply the contrapositive statement: if u{ti) = 
0 then b{ti) = 0. In the applications of this lemma below u will be one of 
the concentrations in a photosynthesis model and the functions a and b are 
obtained by setting the other concentrations to their values in a fixed non¬ 
negative solution. 

The strategy is now to successively obtain restrictions on the position of an 
w-limit point of a positive solution. 

Lemma 2 At an w-limit point of a positive solution of the Pettersson-MA or 
Poolman-MA model on the boundary of S the concentrations of the following 
substances vanish: RuBP, PGA, DPGA, GAP, DHAP, FBP, SBP. 

Proof The proof consists of repeated applications of Lemma I to the solution 
of the dynamical system passing through the w-limit point being considered. 
Gonsider an w-limit point where xqap > 0. Then the evolution equation for 
xpj shows that this quantity is positive at the given point. In the same way the 
concentrations of the following quantities are positive: DHAP, FBP, F6P, E4P, 
X5P, SBP, S7P, R5P, Ru5P, G6P, GIP, ATP, RuBP, PGA, DPGA. This implies 
that the w-limit point is in the interior of S, contrary to the assumptions of the 
lemma. Thus it follows that in fact xgap = 0. The evolution equation for 
xgap then implies that the concentrations of DPGA, DHAP and FBP vanish. 
The evolution equation for xdhap then implies that the concentration of SBP 
vanishes. To proceed further we need to distinguish between the cases where 
the concentration of ATP is non-zero or zero at the given point. In the first case 
we can conclude successively that the concentrations of PGA and RuBP vanish, 
completing the proof. In the second case xp^ = 0 and the evolution equation for 
xpj implies that xpga = 0. It then follows as in the first case that xpuBP = 0. 
Lemma 3 At an w-limit point of a positive solution of the Pettersson-MA or 
Poolman-MA model on the boundary of S the concentrations of X5P, R5P and 
Ru5P vanish. 

Proof The evolution equations for X5P, R5P and RubP show that the con¬ 
centrations of all three substances vanish at an w-limit point of a stationary 
solution if and only if any one of them does. However supposing that none of 
them vanishes leads to a contradiction in the evolution equation for GAP. 
Lemma 4 At an w-limit point on the boundary of S' of a positive solution of 
the Pettersson-MA model with koCA < 5^28 either the concentrations of GIP, 
G6P and F6P vanish or all three are non-vanishing and xatp = xp^ =0. 
Proof In the Pettersson-MA and Poolman-MA models the evolution equations 
for GIP, G6P and F6P show that the concentrations of all three substances 
vanish at an w-limit point of a positive solution if and only if any one of them 
does. In the case of the Pettersson-MA model with ksCA < 5^28 it follows from 
the expression for the time derivative of Li that xatpxgip vanishes at any 
w-limit point. If xgip 7^ 0 at that point then xatp = 0 there. It follows that 
xp^ = 0 at that point 
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Lemma 5 At an w-limit point on the boundary of S' of a positive solution of the 
Poolman-MA model either the concentrations of all three hexoses GIP, G6P 
and F6P vanish or none of them does so. One of the following three cases holds: 

1. xp^ = 0, all hexose concentrations are zero or 

2. xp^ = 0, xatp = 0, all hexose concentrations are non-zero or 

3. xp^ 0, XATP 7 ^ 0, all hexose concentrations are non-zero. 

Proof The first statement is part of Lemma 4. When xgip = 0 at an w-limit 
point on the boundary it follows from the evolution equation for xqip that 
xp^ =0. This is case 1. Otherwise all hexose concentrations are non-zero. In 
that case if =0 the evolution equation for xp^ shows that xatp = 0. On 
the other hand if xatp = 0 the evolution equation for xatp = 0 shows that 

xp^ = 0 . 

Any stationary solution on the boundary satisfies the conditions derived 
above for w-limit points. In the Petterson-MA model a stationary solution with 
non-vanishing hexose concentrations satisfies xp^ = 0 and xatp = 0 without the 
restriction on the parameters occurring in Lemma 4. If, on the other hand, the 
concentrations of the hexoses vanish then the conditions for stationary solutions 
reduce to the condition that either xp^ = 0 or xadp = 0. Given the fact 
that the concentrations of xeap and xs 7 P can be prescribed freely we see that 
there are two three-parameter families of stationary solutions which meet at the 
point where both xp^ and xadp are zero. In the case that the concentrations 
of the hexoses do not vanish we get a three-parameter family of stationary 
solutions satisfying the conditions fc 23 XG 6 P = fc 24 XP 6 P and fesxcip = feexcep- 
Consider now the Poolman-MA model. When xp^ = 0 the two systems agree 
and so the set of stationary solutions is identical. It remains to consider the 
possibility that there are stationary solutions of the Poolman-MA model with 
xpj 7 ^ 0. They would belong to case 3. of Lemma 5 and satisfy the relation 
ksii^CA — Xatp) = ^32- Substituting this into the equation for xp^ gives a 
contradiction and so stationary solutions of this type do not exist. 

4 Linearization about the oj-limit points 

In Lemma 2 and Lemma 3 a set was identified, call it Z, where any w-limit 
point of a solution of the Petterson-MA model or the Poolman-MA model must 
lie. Consider now the linearization of the full model about a point of Z. The 
linearized quantity for a given substance is denoted by a j/ with the name of 
that substrate as a subscript. The aim is to find a block upper triangular form 
for the linearization, so as to obtain information about its eigenvalues. The 
linearized equation for FBP depends only on that substance and so can be split 
off, contributing a negative eigenvalue. At a general point of Z it is difficult to 
proceed further with this strategy. The simplest points about which to linearize 
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are the points zq where all carbohydrate concentrations and xp^ are zero. Denote 
the concentration of ATP at the point zq by a. The linearization is 


dyp, 


■mBP 


dt 

dypGA 

dt 

dyppGA 

dt 

dyATP 

dt 


= —kiypuBP + k22ayRu5P, 

= ‘^^kiypuBP + fc2(cA — a)yDPGA — {k^a + k2s)yPGA, 
= -[fe(cA - a) + k^jyppGA + ksaypGA, 


= fe(cA — a)yDPGA — k^aypGA — k22ayRubP 


-k27ayGiP + k3i{cA - a)yp^, 

dyGAP 


= k^yppGA + keypHAP — kjypAP — k29yGAP, 
= k^yppGA + ki^ysBP + 2fc27a2/GiP + k2?,ypGA 


dt 

dyPi 
dt 

+k29yGAP + ksoypHAP — ksi {ca — a)yPi — k32yPi 

dypHA P 

dt 

dypGP 

dt 

dysiP 
dt 

dyx5P 
dt 

dysBP 
dt 

dysTP 
dt 

dypbP 
dt 

dypubP 

dt 


= —keypHAP + kjyGAP + ki^yspp — ksoypHAP, 
= fesyCeP — fe4yP6P) 

= kisysBP, 

= feo2/fl«5P — k2iyx5P, 

= —{ki3 + ki5)ysBP, 

= kizysBP, 

= kisypudP — ki9yR9P, 

= —kisypubP + ki9yR5p — fc2oyp«5P 


+ k2iyX5P — k22ayRu5P, 
dyG&p 


dt 

dyGiP 

dt 


= —k23yG6P + k24yF6P + fesZ/GlP ” k2eyGep, 
= -fesl/GlP + k26yG6P — k27ayGlP + k32yPi 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 


Here the quantity yppp has been set to zero and its evolution equations omitted 
since it plays no role in what follows. 

Lemma 6 Let L be the linearization of the right hand side of the equations of 
Pettersson-MA model at a point of the form zq. Then all eigenvalues of L have 
non-positive real part and at most four have zero real part. The multiplicity of 
the eigenvalue zero is four when xatp = ca and three otherwise. 

Proof The equation for SBP is decoupled from the others and contributes a 
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negative eigenvalue. The rows and columns corresponding to this variable can 
be discarded. The variables E4P and S7P can then be treated in the same way, 
both contributing zero eigenvalues. Consider next the evolution equations for 
the variables ycip^ VG&p and ypGP- They form a closed system and define a 
submatrix which can be analysed on its own. The trace and determinant of this 
block are negative. The characteristic polynomial can easily be computed and 
the Routh-Hurwitz criterion [3] shows that all eigenvalues have negative real 
parts. The equations for yx 5 P, ypubP and ymp can be treated in exactly the 
same way. Once this has been done the quantity ypuBP can be handled in the 
same way as the variables for FBP and SBP. The quantities ypcA and yoPGA 
define a 2 by 2 block which has negative trace and positive determinant. It 
contributes eigenvalues with negative real parts. Next yGAP and yoHAP can 
be treated together. They also contribute a 2 by 2 matrix with negative trace 
and positive determinant and therefore two eigenvalues with negative real parts. 
The equation for yp. is now decoupled and gives the eigenvalue —k^i^CA — a), 
where a is the concentration of ATP at the given point. Finally yp. contributes 
a zero eigenvalue. This completes the proof. 

Lemma 7 Let L be the linearization of the right hand side of the equations 
of the Poolman-MA model at the point zg with xadp = 0. For ^32 sufficiently 
small there is one eigenvalue with positive real part, three zero eigenvalues and 
all the other eigenvalues have negative real part 

Proof The methods in the proof of Lemma 6 can be used in a very similar way 
to eliminate all the variables except yGiP, yGQP, yp&p^ Vatp and yp. from con¬ 
sideration in the case ^32 yf 0. In the remaining system of five equations yATP 
does not occur on the right hand side and so can also be eliminated, producing 
a further zero eigenvalue. To complete the argument the eigenvalues of a four 
by four matrix must be examined. They depend continuously on the parameter 
fc32. When kg 2 = 0 three of the eigenvalues have negative real parts according 
to Lemma 6. These eigenvalues retain this property for sufficiently small values 
of kg 2 - In this regime the sign of the remaining eigenvalue, which is real, is 
the opposite of that of the determinant. The determinant is — ^24^26^27^32 cA) 
which is negative. 

Lemma 6 and Lemma 7 can be used in combination with the reduction 
theorem (see [7], Theorem 5.4) to obtain results on positive solutions of the 
nonlinear equations which start close to a point of the from zg. Consider first 
the case of the Petterson-MA model with xatp = ca- The centre manifold 
is four-dimensional and is given by the vanishing of all variables except xpiP, 
xs 7 Py Xatp and xp^. The restriction of the system to that manifold is the 
product of a trivial system for xpip and xs 7 P with a two-dimensional system 
which is easily analysed. The conclusion is that solutions which start close 
to Zg have the property that all concentrations converge to limits as t -A 00. 
Generically the limits of xeap and xs 7 p are strictly positive and precisely one 
of the limiting concentrations xadp and xp^ is positive. If instead xatp < ca 
at Zg then the centre manifold is three-dimensional. Again all concentrations 
tend to limits, with the limit of xp. being zero. 

In the case of the Poolman-MA model and the point of the form zg with 
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xadp = 0 there is a one-dimensional unstable manifold and generic solutions 
starting close to Zq do not stay close to zq. Thus zq is stable for the Pettersson- 
MA model and passing to the Poolman-MA model destabilises it. The produc¬ 
tion of sugar from starch in the Poolman-MA model means that there is no 
direct analogue of the argument which was used to show that the function Li 
is a Lyapunov function for the Pettersson-MA model for certain values of the 
parameters. Instead it is possible to identify a parameter regime where —Li is 
a Lyapunov function. 

The total amount of phosphate can be written as the sum of the amount 
contained in the adenosine phosphates, the amount of inorganic phosphate and 
the rest cp = ca + Pr + xp^. For a positive solution the constant ci = cp — ca is 
positive. The quantity Pr can be bounded in terms of Li, counting the carbon 
atoms, with the result that Pr < lOLi. Hence xp^ > ci — lOLi. It follows that 

> 6k32Xp, - C 2 L 1 ( 68 ) 

where C 2 is a constant which depends only on the reaction constants and the 
total amount of adenosine phosphates. Hence 

> 6 ^ 32(01 — lOLi) — C 2 L 1 = 6 fc 32 [ci — (10 -I- 02 /( 6 ^ 32 ))^!]. (69) 

This quantity is positive when ci — (10 -I- 02 /( 6 ^ 32 ) > 0. Let m = 

Then whenever Li is less than m its time derivative is positive. It follows that 
Li is eventually at least m. In particular, it cannot tend to zero for t -A 00 . 
The latter fact follows from the computations of the previous section but here 
an additional quantitative lower bound is obtained for the concentrations of the 
sugar phosphates at late times. It can be concluded that 


liminf(4a:p4P -I- 5xx5P + Qxaip + Qxcgp + bxpgp) > (70) 

t—>oo 

The main results of this section will be summed up in a theorem: 
Theorem 1 There is a non-empty open set of positive initial data for which the 
corresponding solutions of the Pettersson-MA model converge to the boundary 
of S' as t —>■ 00 . They include all data for which the concentrations of all carbo¬ 
hydrate phosphates and inorganic phosphate are bounded by a sufficiently small 
constant e. In the latter case the concentrations of all carbohydrate phosphates 
except E4P and S7P tend to zero while xeap and xrtp remain small as t —> 00 . 
On the other hand, the solutions of the Poolman-MA model corresponding to 
these data do not have the property that the concentrations of all carbohydrate 
phosphates remain smaller than a constant d > 0 for a suitable choice of e. 


5 The Poolman model 

The original model of Poolman uses kinetics which are not mass action. In 
this section we investigate which of the results for the Poolman-MA model have 
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analogues in this case. The total amount of phosphate is conserved and this 
can be derived from the evolution equations independent of the kinetics given 
above together with the fact that {d/dt){xADp) = —{d/dt){xATp)- Lemma 1 
can be applied to this model in exactly the same way as it was applied to the 
Poolman-MA model. For all that is important in that context are the signs 
(positive, negative or zero) of the reaction rates when certain concentrations 
are zero or non-zero. These relations are not changed when the mass action 
kinetics are replaced by the more complicated kinetics in the original model. 
The modulation occurs because the reaction rates are multiplied by positive 
factors depending on the concentrations of other substances. It can be further 
concluded that the analogues of Lemma 2 and Lemma 3 hold for the Poolman 
system. The time derivative of Li is given by 

d 15 

— (5Li) = -V2 - -jVpGA - iVGAP - iVDHAP “ 6 Ust + 6 ^ 17 . (71) 

The manifestly positive term on the right hand side of this equation is given 
by 6 ni 7 = 6 k 32 Xp^. In order to show that the estimate obtained in the case of 
the Poolman-MA model extends to the Poolman model it is enough to obtain 
lower bounds for all the other terms on the right hand side of the equations by 
negative multiples of Li. The quantity V 2 is a combination of two mass action 
terms and there is nothing new to do. The remaining terms admit the desired 
type of bound. 

6 Conclusions and outlook 

In this paper information has been obtained about the dynamics of some models 
for the Calvin cycle of photosynthesis. It was shown for the Pettersson-MA 
model that if the reaction constant determining the rate of conversion of PGA to 
DPGA is small enough compared to that determining the rate of export of PGA 
from the chloroplast then the concentration of some sugar phosphate must attain 
arbitarily small values at late times. In particular for these parameter values 
the model admits no positive stationary solutions. This is a manifestation of 
the phenomenon of overload breakdown. It was further shown that under these 
circumstances all sugar phosphate concentrations except those of E4P, S7P, 
GIP, G 6 P and F 6 P tend to zero for t —>■ oo. For initial data where the initial 
concentrations of xp^, xeap and xstp are sufficiently small all concentrations 
tend to limits at late times and for generic data of this type the limits of the 
concentrations of E4P and S7P are not zero. In contrast it was shown that for 
the Poolman-MA model for generic data with xeap, xsyp and xatp small 
it is not the case that these quantities stay small. Thus, in accordance with 
the original motivation, the regime of overload breakdown is destabilized by the 
additional term introduced by Poolman. This instability is also present in the 
original Poolman model. 

It would be possible to introduce a hybrid model with the network of Pool- 
man and the kinetics as in the Pettersson model. It might be possible to obtain 
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the Pettersson model as a singular limit of this hybrid model but the details 
of how to do this in a mathematically rigorous way remain to be worked out. 
There are a number of other interesting open questions concerning these models. 
Does the Poolman model (or the Poolman-MA model) have a positive station¬ 
ary solution for some values of the parameters? Numerical simulations in m 
indicate that there can be two stable positive stationary solutions which exhibit 
hysteresis when the parameter Vie, which here represents the intensity of light, 
is varied. Compare also the discussion in m- Does the Pettersson model (or 
the Pettersson-MA model) have a positive stationary solution for some values of 
the parameters? Numerical simulations in [5] indicate that this is the case and 
that there can be one stable and one unstable solution. There are, however, no 
mathematical proofs that these features occur. To obtain a positive answer to 
one of these questions concerning the existence of stationary solutions it would 
suffice to show that in the case being considered positive solutions can have no 
w-limit points on the boundary of S and for that it might sufhce to linearize 
about stationary points on the boundary more general than those handled in 
this paper. Beyond the question of the existence of positive stationary solu¬ 
tions there are the questions of their multiplicity and stability. The analogous 
questions for the simple model including ATP considered in and [12] are 
partly open. There are also many other models for the Calvin cycle in the lit¬ 
erature and it would be desirable to do a careful mathematical study of their 
solutions. This is particularly relevant since several errors in the literature have 
been discovered which had gone uncorrected for many years i, m- 
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